BI-TP 2008/12, BNL-NT-08/21 , TIFR/TH/23 



The Spatial String Tension and Dimensional Reduction in QCD 



OO 

o 
o 



M. Cheng^, S. Datta'', J. van der Heide'^, K. Huebner'^, F. Karsch'^'^, O. Kaczmarek'^, 
E. Laermann^ J. Liddle^ R. D. Mawhinney*^, C. Miao^ P. Petreczky'^'^ 
K. Petroyf, C. Schmidf^, W. Soeldner<^ and T. Umeda^ 
'^Physics Department, Columbia University, 
New York, NY 10027, USA 
^Department of Theoretical Physics, 
Tata Institute of Fundamental Research, 
Homi Bhahha Road, Mumbai 400005, India 
'^Fakultdt fiir Physik, Universitdt Bielefeld, 
D-33615 Bielefeld, Germany 
"^Physics Department, Brookhaven National Laboratory, 
Upton, NY 11973, USA 
"RIKEN-BNL Research Center, 
Brookhaven National Laboratory, Upton, NY 11973, USA 

^ DESY, Platanenallee 6, Zeuthen, 
Germany ^Graduate School of Pure and Applied Sciences, 
University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan 



We calculate the spatial string tension in (2+1) flavor QCD with physical strange quark mass 
and almost physical light quark masses using lattices with temporal extent A'^^ = 4, 6 and 8. We 
compare our results on the spatial string tension with predictions of dimensionally reduced QCD. 
This suggests that also in the presence of light dynamical quarks dimensional reduction works well 
down to temperatures l.STc. 



Oh: 



> 

cn 
o 

OO 

o 



PACS numbers: 11.15.Ha, 12.38.Gc 



I. INTRODUCTION 

At high temperature, strongly interacting matter un- 
dergoes a transition from a confining, chiral symmetry 
broken phase to a phase which is chirally symmetric and 
deconfining in nature. For pure gluonic matter as well as 
for QCD with massless quarks this transition is known 
to be a genuine phase transition, while for QCD with 
physical quark masses it is most likely only a smooth 
(though rapid) crossover jH, d, In both cases, ther- 
modynamic observables like energy and entropy density 
show a sharp increase at the transition point, but ap- 
proach the expected behavior of an asymptotically free 
quark-gluon gas only at very high temperatures. 

The phase just above Tc is more complicated than 
a weakly interacting gas of quarks and gluons. Non- 
perturbative features like the (Debye) screening of elec- 
tric modes as well as the possible appearance of a mass 
gap in the magnetic sector give rise to well-known prob- 
lems that make a straightforward perturbative analysis of 
the high temperature phase of QCD fail. While infrared 
divergences arising at a momentum scale p ~ can 
still be taken into account by resummation [1, H, [a, S Hi j 
the softer modes at scale p ~ g^T reflect genuine non- 
perturbativephysics and cannot be handled by any re- 
summation 

The non-perturbative physics arising from the mag- 
netic sector of QCD shows up at distances that are large 
on the scale of the inverse temperature, R ^ l/{g'^T). 
It finds its most prominent reflection in the survival of 



an area law behavior for space-like Wilson loops, i.e. 
confinement of magnetic modes. The spatial string ten- 
sion, as, extracted from the area law behavior of spatial 
Wilson loops, has been studied early on in lattice cal- 
culations performed in the pure gauge sector of QCD 
[Tol [T]| . Detailed studies of the temperature dependence 
of as in SU(2) and SU(3) gauge theories [11 [H, [H 
showed the expected dependence on the magnetic scale 
y/a^ ~ g^{T)T and established a close relation to the 
string tension of a 3-dimensional gauge theory. A recent 
detailed comparison of ^/a^, calculated in lattice regular- 
ized QCD, with calculations performed in the framework 
of dimensionally reduced QCD [15] was extremely suc- 
cessful in describing the temperature dependence of y^oj 
quantitatively even at temperatures close to the transi- 
tion temperature. 

Dimensional reduction makes very specific predictions 
for modifications induced in the magnetic sector of QCD 
through the presence of fermion degrees of freedom, 
quarks. In fact, modifications are expected to be small 
and the dominant structure of e.g. the spatial string ten- 
sion is still expected to be controlled by a purely gluonic, 
3-dimensional SU(3) gauge theory. A possible source for 
failure has, however, been discussed in the past (l6| : dy- 
namical quarks, more precisely quark anti-quark pairs, 
can give rise to additional light bosonic quasi-particle 
modes that are not described by a 3-dimensional gauge 
theory. This makes it interesting to check to what ex- 
tent dimensional reduction also is capable of describing 
physics in the magnetic sector of QCD in the presence of 
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light dynamical quarks. 

In this paper we calculate the spatial string tension 
in QCD with light dynamical quarks. For the first time 
we cover in such an analysis a large temperature range 
and perform calculations at three values of the lattice 
cut-off. We will present a detailed comparison of the 
lattice calculations with predictions based on dimensional 
reduction. Some pre liminary results of this work have 
been presented in |17| . 

This paper is organized as follows. In the next section 
we give a brief discussion of the concepts of dimensional 
reduction as far as they are relevant for this work and 
the analysis of the spatial string tension. In Section III 
we present details of our calculation of spatial pseudo- 
potentials and in Section IV we discuss the determination 
of the spatial string from them. A comparison of these 
results with dimensional reduction is given in Section V. 
Finally we give our conclusions in Section VI. 



II. DIMENSIONAL REDUCTION 

The different non-perturbative length scales of 
0{l/gT) and 0{l/g'^T), respectively, that one encoun- 
ters in the perturbative analysis of QCD at high tem- 
perature are convenient for the reformulation of QCD in 
terms of a hierarchy of effective field theories that clearly 
separate the physics at these different scales 
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allows the disentanglement of sectors that can be treated 
in perturbation theory from sectors that are only acces- 
sible to non-perturbative techniques. 

Since the partition function of an equilibrium field the- 
ory can be written as an Euclidean field theory with the 
time direction being a torus of extent 1/T, the temporal 
modes of fields are proportional to the Matsubara fre- 
quencies, 27rnT, where n is an integer for bosonic and a 
half-integer for fermionic fields. At large temperatures, 
therefore, the non-static bosonic fields and the fermionic 
fields become very heavy and can be integrated out to get 
an effective theory for distance scales^ RT ^ 1 which in- 
volves the static modes only [l^. This gives an effective 
3-dimensional theory^ , 



SE = j d?x\ Tr Fij (x) Fij (x) + Tr [Dj, Ao(x)]^ 
+ m2TrAo(x)2 + AA(TrAo(x)2)2 , 



(1) 



with the static mode of the temporal component of the 
gauge field, A{), appearing as an adjoint scalar field. The 
parameters of this effective theory can be calculated in 
perturbation theory to any order [l^, [23|- At leading 



^ Distances are naturally measured in units of 1/T, and RT ~ 1 
is a natural unit since the average spatial separation of partons 
in the high temperature phase is ~ 1/T. 

^ In this normalization convention the spatial gauge fields have 
dimension one, while the electric field Aq has a dimension 1/2. 



order one has g% = g'^T, = (9 - nf)g*T / {24Tr'^) and 
rnjj = g^T^ {1 + nf/6) for Uf massless quark flavors. 
The deconfined phase of high temperature QCD corre- 
sponds to the symmetric (confining) phase of the above 
3d adjoint Higgs theory (see discussion in [13, HH). The 
confining nature of this theory is the origin of the in- 
frared problem of high temperature QCD for momenta 
p ^ g^T . At very high temperatures also the Debye mass 
will be large, ttid ^ g^T, and the heavy field can thus 
also be integrated out. Therefore, for physics of distance 
scales R » ^/gT, the fields can be integrated out 
which leads to a pure three-dimensional gauge theory. 
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TrFy(x)Fy(x) 
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The gauge coupling g^ can be calculated in terms of gs 
and niD- This has been done to 2-loop accuracy [l^. At 
leading order g^ — g'E- The hierarchy of effective theo- 
ries is referred to as the dimensional reduction scheme. 
The use of the dimensional reduction scheme in combina- 
tion with lattice calculations performed within the effec- 
tive 3-dimensional theory allows the treatment of the in- 
frared problem of perturbative QCD at high temperature 
and the consistent extension of the weak coupling expan- 
sion of thermodynamic quantities to higher orders, for 
instance the 0{g^) contribution to the thermodynamic 
potential ( pressure ) of QCD p3l |. 

In pure SU(N) gauge theory, the underlying nonpertur- 
bative structure of the high temperature theory can di- 
rectly be seen in calculations of spatial Wilson loops. The 
static qq energy in the finite temperature theory is given 
by Wilson loops in space-time planes, W{R,t), which 
are expected to satisfy a perimeter law above Tc, thus 
signaling deconfinement. In contrast, entirely space- like 
Wilson loops, which are restricted to fixed time hyper- 
planes, show an area law, as they can be mapped to the 
qq energy in the confining 3-dimensional theory [l^, Ei 
as indicated above. Detailed studies of spatial Wilson 
loops in SU(2) and SU(3) gauge theories have corrobo- 
rated the dimensional reduction picture [H, [H, . The 
spatial string tension Ug , the coefficient of the area term. 
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has been found to be in very good quantitative agreement 
with that of the dimensionally reduced theory. Since g^ 
sets the mass scale in the 3-dimensional gauge theory, one 
has y^CTs — cg^ where c is a numerical factor. Dimensional 
reduction therefore predicts 



/^s=cg\T)T 



(4) 



to leading order, where g^{T) is the temperature- 
dependent running coupling of the four dimensional the- 
ory and c is a constant that can be determined through 
a non-perturbative (Monte Carlo) calculation of Wilson 
loops in a 3-dimensional SU(3) gauge theory [ll,|23|. Al- 
ready this tree level matching was found to be valid to 
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a good accuracy at temperatures as low as 2Tc [l3|, . 
This is in accordance with studies of other correlation 
functions in SU(2) and SU(3) gauge theories, which also 
show dimensional reduction to be a valid approximation 
already at temperatures T ~ 2Tc. It is worth mentioning 
that the perturbative separation of electric and magnetic 
scales, TfLjj ^ gT > g^T ^ is not valid at these tempera- 
tures, where g{T) > 1 and where the electric gluons thus 
are the lightest degrees of freedom [H, [H, [23| ■ Strictly 
speaking, the step leading from Eq. ([T]) to Eq. ^ thus 
should not be taken. However, the Higgs sector does not 
seem to have a significant effect on the correlation func- 
tions of purely gluonic 3d gauge fields [1^, [27| . 

The situation may be somewhat different in QCD. 
While the quarks acquire an effective mass of 0{7tT), 
they make a qualitative difference in the long distance 
physics: since they provide a source for screening of fun- 
damental charge, one expects that at sufficiently large 
distances RT ^ Stt/^^ (spatial) string breaking takes 
place, which cannot be described by the effective theo- 
ries or 5*3, respectively, as those do not include funda- 
mental charges. Of course, the distance (in units of 1/T) 
one needs to reach to observe this deviation increases 
with temperature, so it remains an important question 
to ask how well the (spatial) string tension of the four 
dimensional theory, extracted at large distances but be- 
fore string breaking sets in, can be explained by the di- 
mensionally reduced theory. In principle it is possible 
to write down the theory which includes also the lowest 
fermion Matsubara frequency 28]. This version of the 
dimensionally reduced theory would be able to describe 
the string breaking as well. At distances smaller than 
the string breaking scale, however, and given the weak 
dependence of the 3d gauge fields on the dynamics of 
the adjoint scalar fields discussed above we would expect 
that the presence of fermionic modes with mass of ©(ttT) 
will have very little influence on Wilson loops, and the 
spatial string tension extracted from them. On the other 
hand, from a comparison of screening masses in a theory 
with Hf — 4 degenerate light flavors, Ref. [31 argued 
that the presence of dynamical quarks could impede di- 
mensional reduction. Thus, the question of an effective 
description of spatial Wilson loops in terms of a dimen- 
sionally reduced theory at not-too-large temperatures is 
not settled and needs to be treated more extensively for 
QCD. 



III. SPATIAL POTENTIAL AT FINITE 
TEMPERATURE 

We have carried out a study of the spatial Wilson loops 
for QCD in the temperature interval [0.9 Tc,4 Tc]. In 
this study we fixed the strange quark mass (m^) to its 
physical value, and the u, d quarks are taken to be de- 
generate with mass O.lm^. These light quark masses are 
about a factor two heavier than physical u, d quarks, and 
give a pion mass of about 220 MeV. The gauge configu- 



rations we analyze here were generated for the study of 
the equation of state of hot strongly interacting matter 
[29| . The gauge action uses Wilson plaquette and rect- 
angle terms. For the fermion action the p4fat3 staggered 
discretization [3Q] was used, which adds a bended three- 
link term to the standard staggered action, and the 1-link 
gauge connections are smeared with the 3-link staples. In 
our analysis we have used lattices with temporal extent 
iV^ = 4, 6 and 8 and aspect ratio N^/Nr > 4. The lat- 
tice spacing has been set using the Sommer parameter rg 
[3H . When quoting the results in physical units the value 
To = 0.469 fm has been used [11]. Given the above tem- 
poral lattice sizes and the temperature interval covered 
in our numerical calculations, the range of lattice spac- 
ings varied from 0.3fm to O.OSfm. Further details on the 
generation of gauge configurations and the action used in 
these calculations can be found in Ref. [29j . 

To interpret the spatial Wilson loops as the expo- 
nential of the static qq potential in the effective 3- 
dimensional theory, we use a transfer matrix formalism 
in the z direction: the loops are constructed as W{R, Z), 
where R is the distance between the static quark and 
anti-quark in the xy-plane. The loops in a timeslice were 
constructed using the Bresenham algorithm ^33| and they 
were averaged over the timeslices. The static qq poten- 
tial, for a transfer matrix acting in the z direction, is 
defined as fSJl 



aVs{R) = - lim In 

Z^oo 



W{R,Z+\) 
W{R,Z) 



(5) 



To obtain a better signal-to-noise ratio it is impor- 
tant to be able to reach the asymptotic limit already for 
moderate values of Z. To achieve this, the gauge fields 
connecting quark and anti-quark were smeared using 3- 
link staples perpendicular to the direction of the gauge 
field. The smearing procedure has been applied itera- 
tively. We found that for small values of the gauge cou- 
pling P {(3 < 3.41) 10 APE smearing steps with a smear- 
ing coefficient of 0.4 were sufficient to give the asymp- 
totic value already at Z = 3. For larger (3 we use 30 to 
60 smearing steps. 

The Wilson loop ratios in Eq. ^ were performed 
for different and a plateau is reached by Z=3. To 
avoid any remaining contamination from higher states, 
the logarithm of the ratio of Wilson loops was fitted to 
the ansatz —aVg + b ■ exp{—cZ). The results from the 
direct ratio at Z — 3 generally agree well with the value 
obtained from fits. We used in our analysis estimates for 
the string tension coming from Eq. ([5]) with Z = 3 as well 
as results obtained from an extrapolation to = cxd. Dif- 
ferences in these results have been treated as a systematic 
error in our determination of as- 

At each value of the temperature we have generated 
several thousand (up to 33,000) gauge field configura- 
tions using a Rational Hybrid Monte-Carlo (RHMC) al- 
gorithm. The total number of gauge configurations gen- 
erated at each gauge couphng can be found in Ref. [29| . 
The calculation of spatial Wilson loops has been per- 
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formed on gauge configurations separated by 10 RHMC 
trajectories. To reduce autocorrelation, the measure- 
ments were bunched into blocks and errors on aVs{R) 
were calculated using Jackknife samples. By varying the 
block-size of the Jackknife samples we checked that the 
calculated errors do not change significantly. 
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FIG. 1; The spatial potential Vs{R,T) calculated on A^r = 4 
(top) and A''t — 6 (bottom) lattices. The reduced temperature 
T/Tc has been calculated using the values of ro/a determined 
in [291 as well as the values roT^ = 0.484 for AT^ = 4 3 and 
roTc = 0.466 for iVr = 6 [l!]. The line is the fit to zero 
temperature potential calculated at /3 = 3.63. 

Let us first discuss the most prominent features of the 
spatial potential at finite temperatures and compare it 
to the usual zero temperature potential. Note that for 
the zero temperature lattices, of course, timelike and spa- 
tial Wilson loops are identical. In Figure [1] we show the 
spatial potential calculated for Nt = 4 and Nt = 6 at 
a few selected values of the temperature. We have sub- 
tracted the renormalization constant determined at zero 
temperature in Ref. |29j to make the comparison between 
different temperatures easier. Figure [T] shows that be- 
low Tc the spatial Wilson loop shows very little change 
with temperature, while above the transition tempera- 



ture a strong temperature dependence is seen. Only at 
very short distances the spatial potential remains, within 
errors, temperature independent also above Tc- For tem- 
peratures close to the transition temperature the spatial 
potential falls below the zero temperature potential. The 
effect is more pronounced for Nr = 4 than for Nt — 6 
which might indicate that this is a cut-off effect. As the 
temperature increases the slope of the spatial potential at 
large distances, i.e. the spatial string tension, clearly in- 
creases in agreement with observations made previously 
in pure gauge theories. This will be discussed in more 
detail in the next section. 

We find no evidence of string breaking occurring in the 
high temperature spatial potentials. This, however, may 
not be too surprising as it is known from studies of the 
heavy quark potential at zero temperature that smeared 
Wilson loops at numerically accessible perimeters are not 
well suited for studying string breaking. Polyakov loop 
correlation functions have been found to perform much 
better in this respect [sHi- Therefore one might want to 
also analyze spatial Polyakov loop correlation functions 
to get more sensitive to string breaking effects in spatial 
potentials. 



IV. SPATIAL STRING TENSION 

In order to extract the spatial string tension from the 
large distance behavior of the spatial pseudo-potentials 
we have to rely on fits. However, the choice of an appro- 
priate fit Ansatz is less obvious than at zero temperature. 
At zero and low temperatures the (pseudo-) potential is 
usually well described by 



Vs{R) = -- + co 



(6) 



The l/R term, in 3 + 1 dimensions, may arise from 
two contributions, from Coulombic perturbative gluon 
exchange at small distances, approximately <0.2 fm, and 
from string fiuctuations relevant at distances >0.3 fm. 
The later is often referred to as the Liischer term. In 
fact, in the proportionality factor of the string fluctuation 
term, a^, is a universal constant. 



{D - 2)tt 
24 



(7) 



that only depends on the space-time dimensionality D 
[36| . As dimensional reduction manifests itself as temper- 
ature increases one should expect that the contribution 
arising from string fiuctuations should gradually turn 
from 7r/12 to 7r/24. In fact, finite temperature correc- 
tions to the string fluctuation term have been calculated 
previously [s^l and have been found to be of relevance for 
the analysis of the temperature dependence of the (con- 
ventional) heavy quark potential [38|. At the same time 
and for the same reason, the Coulombic term due to per- 
turbative gluon exchange is expected to change from a 
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1/R behavior appropriate in 4 dimensions to a logarith- 
mic R dependence in an effectively 3 dimensional theory. 
It is this gradual change in the short and intermediate 
distance part of the pseudo-potential that is responsible 
for ambiguities in the fit Ansatz. We will discuss in the 
following our strategy to deal with this ambiguity. 

From our data, we could not clearly disentangle a 
log(i?) term from a 1/R behavior. In fact, addition of 
the log(i?) term to the form of Eq. ([6]) only makes the fit 
very noisy. Similarly we could not restrict the fit inter- 
vals to such large distances where 1 /R terms of whatever 
origin, Coulombic or string fluctuations, could be safely 
neglected. We therefore have fitted our data to Eq.® 
with 3 free parameters. We note that this is not too se- 
vere a restriction as we are not interested in a detailed 
description of the i?-dependence of the pseudo-potentials 
but rather want to get control over its large distance be- 
havior. 

On our finer lattices {N^ — 6 and 8 ) we observe that 
the coefficient of the term, i.e. a, decreases with in- 
creasing temperature and approaches 7r/24 at high tem- 
peratures. This is shown in Figure O The fit results for 
a are somewhat dependent on the fit range. We have 
adopted to choose a minimum separation Rmm between 
0.7 and 1 times tq. At sufficiently large distances the 
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FIG. 3: The string tension in units of ro as function of the 
temperature calculated on A''t = 4, 6 and 8 lattices. For 
better visualization the temperature axis has been scaled by 
To = 2mMeV corresponding to roTo = 0.47619. 



The values of the spatial string tension for different A',- 
are summarized in Table U For A^,- = 6 and 8 the string 
tension has been determined from 3 parameter fits. For 
the A^T- = 4 case we used mostly two parameter fits with 
different values of a. The errors on the spatial string 
tension shown in the table are predominantly systematic 
due to the choice of the fit form and fit interval. The 
temperature dependence of the spatial string tension is 
shown in Figure [31 

For temperatures T < 1.2Tc the spatial string ten- 
sion is close to the zero temperature string tension and 
is linearly rising with temperature for T > l.bTc This 
is expected if dimensional reduction holds at these tem- 
perature values. We will discuss this in more detail in 
the next section. It is possible that also the rapid drop 
in the coefficient a of the Coulomb term for T > l.STj. 
is related to the onset of the 3-dimensional physics re- 
ficcted in the linear rise of the spatial string tension in 
this temperature region. 



FIG. 2: Fitted values for a compared with the coefficients of 
the Liischer term for 3 and 4 dimensions, Nt — 6, 8. The 
temperature has been scaled by To = 20071^61^, where roTo = 
0.47619. 

Coulombic term will not be seen, with only string fluctu- 
ations contributing. Fixing the value of a to the coeffi- 
cient of the Liischer term, ul, may describe the potential 
at large distances. This Ansatz was fitted to the poten- 
tial at distances greater than tq. For small temperatures, 
T < 1.5Tc, with a fixed to 7r/12 the values of a coming 
from these 2 parameter fits are systematically larger than 
those coming from the 3 parameter fits. At high temper- 
atures, however, the spatial string tension obtained from 
two parameter fits with a = 7r/24 and three parameter 
fits agree within statistical errors. 



V. COMPARISON WITH THE PREDICTION 
OF DIMENSIONALLY REDUCED THEORY 

As has been discussed in Section II we expect that at 
high temperatures the spatial string tension should be 
given by y/a^ — c • g%{T). At temperatures several or- 
ders of magnitude larger than the transition temperature, 
when the static electric field Aq can be integrated out, 
the effective theory is just a pure SU(3) gauge theory. 
The proportionality coefficient c is just a constant and 
has been determined to be [l^l c = 0.5530(10), corrobo- 
rating an earlier value of 0.554(4) [l3|. In the interesting 
temperature range of a few times Tc, however, c will de- 
pend on the mass tod and coupling \a of the Aq field and 
thus on the temperature. In the light of the approximate 
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TABLE I: The spatial string tension in lattice units extracted from Coulomb plus linear fits on = 4, 6 and 8 lattices (see 
text for the details). The temperature scale, roT, is obtained from the values of ro/a calculated in Ref. [29l ]. 



decoupling of the 3d scalar and gauge fields we expect 
that the dependence of the coefficient c on these param- 
eters should be weak and its value should be close to the 
pure gauge value given above. Indeed, the calculations 
of the string tension in 3-dimensional adjoint Higgs mod- 
els show only a weak dependence on the parameters of 
the scalar sector and its value is close to the pure gauge 
value [l^]- Unfortunately, the statistical accuracy of the 
spatial string tension calculated in the 3d adjoint Higgs 
model is significantly lower than in the pure gauge case 
and no continuum extrapolation has been performed. For 
the relevant case of three quark flavors and temperatures 
of about 2Tc the value of c is about 2% lower than the 
pure gauge value at fixed lattice spacing corresponding to 
Q/{g^a) — 21. The analysis of Ref. [l^l also suggests that 
a possible temperature dependence of c is less than 5%. 
Therefore we use an averaged value c = 0.54(1) in our 
analysis, where the indicated uncertainty also includes 
possible temperature dependence in it. 

The gauge coupling of the effective 3d theory has been 
calculated to 2-loop accuracy [lB| 



a\n) + 9\n) (2^0 In 

g%f.) (25iln(A) + 



02 + 



26o In ( 



f ^J■ \ 8n/ In 2 — 3 



\t:TJ 



487r2 



«3 (8) 



Here g{fi) is the QCD couphng in MS scheme, 6o 
(11 - 2n//3)/(47r)2 and h = (102 - 38n//3)/(47r)4 g 



the coefficients of the universal 2-loop beta function, and 
p, is the MS renormalization scale. The coefficients 02 
and 03 can be found in Ref. p^ . In our case n/ = 3. 
The coupling depends on the renormalization scale 
p, at any fixed order of perturbation theory. Of course, 
the dependence on jj, gets weaker and weaker as we go 
to higher orders of the perturbative expansion. In prac- 
tice, however, we have to deal with the scale dependence 
of the effective coupling. Following Ref. [l^ we fix the 
scale /I, using the principle of minimum sensitivity, i.e. 
we require that the derivative of the 1-loop expression for 
g^{T,p) vanishes at /l = /l*, and vary the scale in the 
interval jl ~ (0.5 — 2)/i*. For 3-flavor QCD we find the 
value /i* ~ 9. IT. To specify g{p,) and thus gE completely 
we need to know the value of A^g^, or more precisely the 
ratio T/ Ajjg. Since the temperature has been set by the 
Sommer parameter rg this means that we need to spec- 
ify AjjgTo. This can be done in principle by calculating 
the ro parameter at several lattice spacings and fitting it 
by the modified 2-loop Ansatz [s^ to determine rQAj^at- 
One then may use lattice perturbation theory to calculate 
A-p g/A^gf. This has been done in SU(3) gauge theory 
|4C1 | as well as in 2-flavor QCD with Wilson fermions [4l| . 
Unfortunately, the perturbative calculations needed for 
this have not been performed for the p4fat3 action. On 
the other hand we can express the temperature in phys- 
ical units using the Sommer parameters tq = 0.469(7)fm 
obtained from quarkonium splitting as input [3^ . Fur- 
thermore, using the same physical input the running cou- 
pling constant has been calculated on the lattice in 2-1-1 
flavor QCD within the so-called V scheme [i^l, giving 
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0.2082(40). 
2!(/i)/(47r) 



This corresponds to 
0.180(+2)(-3) if we 



av{n = 7.5GeV) 
"ms(/^ = 7.5GeV) 

use the 3-loop relation [43| between the coupHng in the V- 
scheme and the MS'-scheme. Using the 2-loop beta func- 
tion we can determine the couphng g{fl) entering Eq. ([8]) 
at the given scale /i* thus specifying the gauge coupling 
of the effective 3d theory g%{T) for different tempera- 
tures. Combining the value of c from above and g%{T) 
we finally get the corresponding prediction for the spa- 
tial string tension in the dimensionally reduced theory. 
In Figure m we show our results for T/ ^Jo\JT) compared 
with the prediction coming from the dimensionally re- 
duced theory as discussed above. This is shown as a 
solid line with its uncertainty shown as a band (dashed 
lines) . There are three sources of uncertainty in the di- 
mensional reduction prediction. The first is the uncer- 
tainty in the value of c. This is the dominant source of 
the error at high temperatures, T > 2Tc. The second 
source of error is the scale dependence oi gE, which is 
the most important one at low temperatures, T < 2Tc. 
Finally there is an error in g{fi) coming from the value 
of ajjg{7.5GeV) ^ 0.180(+2)(-3). This, however, is 
significantly smaller than the previous two in the entire 
temperature range. 

Figure m seems to suggest that dimensional reduction 
works down to temperatures surprisingly close to the 
transition temperature. In order to see whether this pic- 
ture is self-consistent one has to calculate other spatial 
correlation functions. In particular, one has to verify 
that the largest correlation length of operators built from 
quarks is sufficiently small to justify integrating out the 
Matsubara modes of quarks. In Ref. this problem has 
been studied in 4-flavor QCD, using Nt = 4 lattices with 
the standard staggered formulation. The analysis per- 
formed in that paper suggests that the pion correlator 
gives the largest correlation length up to temperatures 
as high as 3Tc. Note, however, that this may be due to 
large discretization errors in the screening masses when 
the standard staggered action is used on = i lattices. 
Indeed, for this action it has been noticed in Ref. 
that the pion screening masses on Nr = 10—12 lattices 
may come close to a value of 2ttT. Recent preliminary 
calculations with improved (p4fat3) staggered fermions 
[45, 46] give large values for this quantity, rriT^/T ~ 4.6 
for T ~ 1.5Tc and m„/T ~ 5 for T ~ 2.0Tc, aheady 
on lattices with Nr = 4 and 6. Thus the value of the 
pion screening mass is larger than the smallest glueball 
screening mass for T > l.bTc which was estimated to be 
~ AT [16i | and larger than the Debye screening mass es- 
timated in This indeed would suggest that in QCD 
dimensional reduction may work down to temperatures 
as low as l.STc also in the case of 2-1-1 flavor. 



VI. CONCLUSIONS 

In this paper we have calculated the spatial string ten- 
sion in 2-1-1 flavor QCD with a physical strange quark 
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FIG. 4: T/^Os{T) calculated on iV^ = 4, 6 and 8 lattices 
compared with the prediction of dimensional reduction indi- 
cated by the line. The uncertainty in the prediction of the 
dimensionally reduced theory is shown by dashed lines. For 
easier visualization the temperature on the lower horizontal 
axis has been rescaled by To = 200 MeV corresponding to 
roTo = 0.47619. We also show the temperature scale in units 
of ro on the upper horizontal axis. 



mass and light quark masses of 0.1ms corresponding to a 
pion mass of about 220 MeV. The spatial string tension 
calculated at different lattice spacings agree reasonably 
well with each other. We have compared the results of 
our calculation with the prediction of dimensionally re- 
duced effective theory and have found remarkably good 
agreement down to temperatures close to the transition 
temperature. This is similar to the observation made in 
SU(3) gauge theory [l^. There are three sources of un- 
certainty when comparing the data on the spatial string 
tension the scale dependence of the 3d gauge coupling, 
the uncertainty in the value of coefficient c and the un- 
certainty in roA-p-g. The uncertainties from the last two 
sources could be reduced by calculating the lattice beta 
function for the p4fat3 action perturbatively and through 
a more precise calculation of the string tension of the 3d 
adjoint Higgs model. 

Let us finally note that the spatial string tension has 
been recently studied also in 2 flavor QCD using Wilson 
fermions and significantly larger quark masses [H, |4^ . In 
Ref. [4^ the spatial string tension has been calculated 
only up to 1.28rc and no temperature dependence has 
been found in this temperature interval. The results of 
Ref. [495 , obtained on lattices with temporal extent Nr — 
4 qualitatively agree with ours, but the drop of the string 
tension close to is significantly larger. It remains to 
be seen to what extent these discrepancies are due to 
larger quark mass, cutoff effects or limited statistics, as 
calculations with Wilson fermions are numerically more 
demanding. 
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